Entanglement Entropy of the Two-Dimensional Heisenberg Antiferromagnet 
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We compute the von Neumann and generalized Renyi entanglement entropies in the ground-state 
of the spin- 1/2 antiferromagnetic Heisenberg model on the square lattice using the modified spin- 
wave theory for finite lattices. The addition of a staggered magnetic field to regularize the Goldstone 
modes associated with symmetry-breaking is shown to be essential for obtaining well-behaved values 
for the entanglement entropy. The von Neumann and Renyi entropies obey an area law with 
additive logarithmic corrections, and are in good quantitative agreement with numerical results 
from valence bond quantum Monte Carlo and density matrix renormalization group calculations. 
We also compute the spin fluctuations and observe a multiplicative logarithmic correction to the 
area law in excellent agreement with quantum Monte Carlo calculations. 
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I. INTRODUCTION 

The scaling of entanglement entropy in the ground- 
states of quantum many-body Hamiltonians has be- 
come a useful tool for characterizing certain universal 
propcrtiesJj For the ground-state wavefunction |^) the 
entanglement entropy is defined as the von Ncumman 
entropy of the reduced density matrix pj2 = Trj^j^') (vl>| 
of subsystem fl obtained by tracing out the degrees of 
freedom in the remainder of the system fl, 



Sipn) -Tr(pj2ln/5o) 
More generally, the a-Rcnyi entropies are defined as 

1 



Saipn) = 



1-a 



(1) 



(2) 



and reduce to the von Neumann entropy in the limit a — > 
1. In one dimension the behavior of the von Neumann 
and Renyi entropies is understood quite well, especially in 
critical systems which, exhibit universal logarithmic scal- 
ing with system size, El but in higher dimensions the von 
Neumann and Renyi entropies are difficult to study. Even 
numerically, determination of the entanglement entropy 
is limited to exact diagonalization for small system sizes 
and the density matrix rcnormaliMtjon group (DMRG) 
for quasi-one-dimensional systems.&cl Recently the Renyi 
entropy 52 was calculated using qusintum Monte Carlo 
(QMC) in the valence bond basis,El but the von Neu- 
mann entropy remained inaccessible in QMC. In dimen- 
sions d > I some general results are available for 
behavior of tha entanglement entropy in free-fermior 
and free-bosonH theories, and in two dimensions universal 
terms p^¥^ also been identified for topologically-ordcred. 
phaseaij'EJ and conformal quantum critical points£3'L3 
Detailed calculations for specific microscopic Hamiltoni- 
ans are few, however. 

In this work we compute the von Neumann and gen- 
eralized Renyi entropies of the spin-1/2 Heisenberg an- 
tiferromagnet on the square lattice using the modified 



spin- wave theory for finite latticesll3Tl3 and show that, 
as expected, thc_.von Neumann and Renyi entropies 
obey an arca-JaMlZI with important additive logarithmic 
corrections.E3li3li3 For the case of the square lattice with 
periodic boundary conditions (PBCs) most of the calcu- 
lations can be done analytically, with the exception of 
finding the eigenvalues of a matrix. For ladder systems 
with open boundary conditions (OBCs) we present a sim- 
ple numerical method that generalizes the procedure for 
the translationally-invariant case. We performed QMC 
and DMRG simulations and show that the spin- wave re- 
sults compare quite favorably with these numerical find- 
ings; in addition, thfre is excellent agreement with recent 
numerical results.HD. We also compute the spin fluctua- 
tions and observe a multiplicative logarithmic correction 
to the area law in excellent agreement with QMC results, 
further supporting the accuracy of the spin- wave theory. 

Although Rcf. P mentions a spin-wave calculation for 
the entanglement entropy of Heisenberg ladders, this is 
the first time a proper treatment of the Goldstone modes, 
which strongly affect the entanglement entropy, has been 
carried out; previous spin- wave calculations have also 
missed the additive logarithmic correction. 



II. MODIFIED SPIN- WAVE THEORY 

We begin with a careful treatment of the spin-S* antifer- 
romagnetic Heisenberg model on a d-dimensional hyper- 
cubic lattice with sublattices A and B, which is described 
by the Hamiltoniann3 



(3) 



with J > 0, where S are the vectors connecting a site 
on one sublattice with its z — 2d nearest neighbors on 
the other sublattice, and (—1)'' = 1 for r £ A and — 1 
for r G B. We assume PBCs. The staggered field h > 
has been included to avoid spurious divergences in fi- 
nite lattices. Physically, the staggered field regularizes 
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the Goldstonc modes associated with the spontaneous 
breaking of the 0(3) spin-rotational symmetry of Hamil- 
tonian (^) for = 0. Indeed, although the spin- wave 
theory is an expansion around the classical Neel state in 
which spins on sublattice A have magnetization S and 
spins on sublattice B have magnetization —S, we will 
proceed to restore the sublattice symmetry by adjusting 
h such that (S^) — 0. The result of this modified spin- 
wave theory, treated as a variational wavefunction for 
Hamiltonian (^) with h = 0, has been shown to give sensi- 
ble and remarkably accurate answers for many quantities 
of interest as-eompared to quantum Monte Carlo (QMC) 
calculations .Ej In this work we show that the spin- wave 
theory also gives accurate answers for the von Neumann 
and Renyi entanglement entropies, as well as the fluctu- 
ations of spin . 

We use the Dyson-Maleev representation to write the 
spin operators in terms of bosonic creation and annihila- 
tion operators 6j, br as 5*+ — {2S~hr)br, S~ = &J, 5*^ ~ 
S-hr for r e A and S+ = -bl{2S-nr), = -K, = 
— {S — hr) for r G i?, with n,. = bib,.. In terms of the 
bosonic operators we have for r £ A 



S{nr 



n-r+s 



— brb 



rOr+S 



bib's) 



'-,bl(bl+s-br)X+s, 



(4) 



so that, after dropping the fourth-order term in Eq. (jj) 
(equivalent to making a single-particle density matrix 
ansatzo), Hamiltonian becomes the linear spin- wave 
Hamiltonian 



LSW 



E^eel + {zJS +h) 



^J2ibrbr+S+h.C.) 
^ r.S 

(5) 

where Ej^cci/N = —zJS'^/2 — hS is the energy of the clas- 
sical Nccl state. The truncation of the bosonic Hamilto- 
nian has the effect of introducing small mixing between 
the physical states with fi^ < 25* and unphysical states 
with hr > 2S, which we will see leads to a slight over- 
estimation of the entanglement entropy due to the larger 
number of degrees of freedom forming the ground-state. 

Although it is conventional to diagonalize Hamilto- 
nian treating the two sublattices separately, we can 
simplify many of the formulas below by treating all 
bosonic operators as equals. Introducing the Fourier 
transform b^ = (1/v^) X^k '*""^k where N = L'^ is 
the total number of sites and fc^ = 27rn^j/L for = 
—L/2 + 1, . . . , L/2 in each direction /i = 1, . . . , d, Hamil- 
tonian (H) becomes 



LSW 



Nccl ■ 



k 



1] 



"-k ^(OkO-k + h.c. 



with rik ~ blJj^ and 



h 

7js 



7k 



ik (5 



(6) 



(7) 



Using the canonical transformation /3k = cosh 0k ^k 



sinh^k&Lk with tanh(26'k 



_ = So H 
2 and Eq = 



finally written as -ffLSW 
^ {zJS /t])^1 - (rnk) 
X^k '^k]/2 (this is not the variational ground-state energy, 



Hamiltonian (|^) is 
Ek'^k^k/^k where 

- -EnooI - [zJSN/rj- 



however, see Ref. 
values are 

1 



1141 ). The single-particle expectation 



(fet&rO = --<5rr'+/(r-r'), {b,b,.) ^ g{v - v') , (8) 



where 



/(r) = TTTT cos(k • r) — , 

9(r) = TTTT cos(k • r) — ■= 
^ ' 2N ^ ^ ' ,/T 



v/1 - (^7k)2 



(9) 



(10) 



In particular, e — {hr) = /(O) — 1/2. We now focus on 
the square lattice, d = 2. Typically, we are interested 
in only the thermodynamic limit N ^ oo for which the 
staggered magnetization order parameter at /i = is 



rriz = lim lim 



-E(- 



{SI 



s- 



(11) 



where Coo — 0.197 after converting Eq. (^) into an inte- 
gral. For finite TV, however, the k = (0, 0) and k = (tt, tt) 
Goldstone modes corresponding to the breaking of 0(3) 
spin-rotational symmetry on each sublattice contribute 
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1 



1 



1 



yJh/2JS 



0{^Jh/2JS) 



(12) 



to /(O), which shows explicitly how e, and hence m^, 
diverges for /i — > unless the limit — oo is taken 
first. This is of course the order of limits required for 
a spontaneously broken symmetry. We can "repair" this 
divergence for finite lattices within the spin- wave approx- 
imation by requiring that =r-0, i-c, by adjusting the 
staggered field h so that e = S':til 



/(O) -\ = s. 



(13) 



This amounts to restoring the broken sublattice symme- 
try by hand. As expected, /i — > as A — > oo, and indeed 
for large A wc have 



1 1 



AJS 2ml ' 



^> 1. 



(14) 



Various observables calculated within this formalism 
have been shown to agree quite well with QMC calcu- 
lations on finite lattices, withrii.i.t the need to extrapolate 
to the thermodynamic limit.oEj As we will see below, 
this property is also essential for observing the correct 
scaling behavior of the entanglement entropy and spin 
fluctuations with system size. 
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III. SPIN FLUCTUATIONS 

Eq. ( |l3|) is the key step to finding well-behaved values 
for the entanglement entropy for finite system sizes. In 
particular, it is not sufficient to use small but arbitrary 
values of h and extrapolate to ft, — 0. To emphasize this 
point, we first consider the spin fluctuations in a region 

n, 



.F(f2) = ((5f,-(5f,)n, 



(15) 



The fluctuations have been shown to exhibit scaling 
behavior similar (though not identical) to the entan- 
glement entropy in several systems, including most 
gapped piodjcls and notably in onc-dimcnsional critical 
systcmsE2l'c3 In free fermion-likc models in any dimen- 
sion they are reia.ted directly to the von Neumann and 
Renyi entropdesJEao While for other models the relation 
is less directE§o and in some cases non-existent ,lZI we 
can nevertheless use J^i^), which can be computed eas- 
ily in QMC calculations, as a first check on whether the 
ground-state obtained in the spin-wave calculation cor- 
rectly captures non-local correlations. 

The spin fluctuations are computed by summing the 
spin-spin correlation function as 



(16) 



r.r'efi 



where 



' {{nrflr') - (fir) {nr')) (17) 



with Crr' = 1 if r, r' arCj^ the same sublattice and — 1 
otherwise, and, as usual,La a factor of 1/3 has been intro- 
duced to compensate for the loss of spin-rotational sym- 
metry, (Sr ■ Sr') = (S^S^i) in the present approach. In 
other words, we consider (Sr • Sr') as the primary quan- 
tity and derive (S^S^,) from it. Using Eq. (||) and Wick's 
theorem, we find {hrfir') — ("■r)(f^r') = —Srr'/4: + f{r ~ 
1"')^ + .9(1" ~ i"')^- Since g{r — r') vanishes when r, r' are 
on the same sublattice while /(r — r') vanishes when r, r' 
are on different sublattices, we obtain 



{s:s:,)-{s:){s^;, 
_ 1 ' 
~ 3 



+ f{v-v'f-g{v-v'f 



(18) 



The factor of 1/3 in Eq. (|17|) requires some further 
justification. The spin-wave theory is highly anisotropic, 
which should not be the case if, for finite lattices at ft, = 0, 
we demand that the spin-rotational symmetry be pre- 
served. For example, we recall that the exact ground- 
state of Hamiltonian (^ commutes with the total spin 
operator in the z-direction S^^^ = , which implies 

that only the diagonal elements of the reduced density 



0.15 



is 




0.05 



FIG. 1: (color online). Ground-state spin fluctuations in the 
spin- 1/2 Heisenberg antiferromagnet on L x L square lattices 
with PBCs. For each L the subsystem is half of the total sys- 
tem, i.e., an I//2 X L region as shown in the inset. Dashed lines 
show fluctuations computed in the spin-wave approximation 
for different L at fixed values of the staggered magnetic field. 
Solid lines show fits to the form J-{L) — aLlnL + bL + c, 
with a ~ 0.023, b ~ 0.065 for the spin-wave results and 
a ~ 0.023, b ~ 0.066 for the QMC resuhs. 



matrix for a single spin are non-zero and for spin-1/2 
take the form 



Pr = ' 2 







(19) 



With the condition that (S^) = 0, we expect ((5^)^) = 
1/4, while apparently {(S^)'^) = ((Sr)^) = S{S + 1) = 
3/4. Thus the proper interpretation is that we first com- 
pute ((Sr)2) and use the expected relation {{S'^)'^) = 
((5-^)2) ^ ((5-^)2) to determine {(S^)^), and similarly for 
more general correlation functions. The factor of 1/3 is 
also consistent with the expectation that spin-wave the- 
ory slightly underestimates fluctuations relative to QMC 
as in Fig. n^. 

QMC calculations for the Heisenberg antiferromagnet 
show that J^(r2) obeys an area law with a multiplicative 
logarithmic correction, but for any fixed value of the stag- 
gered magnetic field the spin-wave calculations suggest a 
strict area law as shown in Fig. |l|, and moreover, extrap- 
olating to ft — > docs not give the desired result. In 
contrast, the fluctuations computed using the prescrip- 
tion (S^) = show excellent agreement with QMC data, 
especially for the fit to the form J^{L) = aL\n L + bL + c. 
We note that a similar multiplicative logarithmic correc- 
tion was also observed for the valence bond entanglement 
entropy, which counts the number of singlet bonds cross- 
ing the bound|ax5z_between the subsystem and remainder 
of the system.EaS 

The close agreement between the spin- wave and QMC 
results gives evidence that the ground-state of Hamilto- 
nian (^ computed in the modified spin- wave theory ade- 
quately accounts for non-local correlations as quantified 
by the spin fiuctuations. 
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FIG. 2: (color online). Spin- wave results for the ground-state 
Renyi entropy S2 in the spin- 1/2 Heisenberg antiferromagnet 
on LxL square lattices with PBCs. For each L the subsystems 
are £ x £ squares as shown in the inset, cf. Ref. m. 



IV. VON NEUMANN AND RENYI 
ENTANGLEMENT ENTROPIES 




L 



FIG. 3: (color online). Ground-state von Neumann entangle- 
ment entropy 5 and Renyi entropies S2, Ss, S4 in the spin- 
1/2 Heisenberg antiferromagnet on L x L square lattices with 
PBCs. The setup is the same as Fig. |l|, i.e., subsystems of 
size L/2 X L. Dashed lines show the von Neumann entropy 
computed in the spin-wave approximation for different L at 
fixed values of the staggered magnetic field. Solid lines show 
fits to the form S{L) = aL + bin L + c, cf. Fig. § 



With the fluctuations showing excellent agreement 
with QMC results in Fig. ^ we now determine the en- 
tanglement entropy within the modified spin-wave the- 
ory. Let Xrr' = {{br + St)(6r' + bl)) /2 and Prr' = 

-{{k-bl)iK,-bl,))/2,or 



Xrr' = /(r-r')+.g(r-r'), 
Prr' = /(r - r') - .g(r - r'). 



(20) 
(21) 



Following Ref. ^ (and the equivalent method in Ref. 31 ), 
the entanglement entropy of a region n can be computed 
from the eigenvalues of the matrix 



'Pr' 



(22) 



r"eo 



where r,r' e ft. It is important to note that Crr' is not 
necessarily symmetric, and a suitable method for finding 
the eigenvalues must be used. The eigenvalues can be 
shown to be real and satisfy > 1/4, however, in accor- 



dance with the uncertaintji 
entropy is then given byE2H 



principle. The entanglement 



Sin) 



E 

9 



-] In ( l^a + - 
2 \ ' 2 



1^1 f 1 



(23) 



In particular, when the subsystem is a single spin, 
Crr' = /(O)^ SO that from the constraint (^3|) we have 
5(1) = (1 + 5')ln(l + S) - S\nS. For spin-1/2 this 
gives 5(1) ~ 0.955, which actually over-estimates the 
exact answer given by the reduced density matrix (19), 
ln2 ~ 0.693. Although particularly extreme for very 



small subsystem sizes, a slight over-estimation is found to 
be a general feature of the von Neumann entropy com- 
puted in the spin-wave theory, and is most likely due 
to entanglement arising from the unphysical degrees of 
freedom that mix with the physical degrees of freedom 
when the original Hamiltonian is truncated. We can also 
compute the Renyi entropies, which are given by 



5,(0) = 



+ 1/2)" - {u, 1/2)"]. (24) 



The second Renyi entropy is particularly simple, 
S2{n) =p-^« ln(2!^g), as is the single-copy entanglement 

cntropy&Ea 5oo(f^) = X]g ln(i^g-l-l/2). For 52 the single- 
site Renyi entropy from the spin-wave theory is exact 
for spin-1/2: ^2(1) = In 2, the same as would be from 
Eq. (^9|); this is special to ^2, however. The Renyi en- 
tropy ^2 provides another important check on the validity 
of the spin- wave approximation in computing the entan- 
glement entropy, as it is not clear a priori that such a 
simple theory should properly capture the entanglement 
content of the ground-state. But 52 was recently com- 
puted in valence bond QMC, and as shown by a com- 
parison of Fig. 1^ with Ref. ^ the spin- wave results are in 
remarkably good agreement with numerical data. Spin- 
wave theory again slightly over-estimates 52 . Fig. ^ also 
presents 52 for much larger total system sizes than was 
possible in Ref. ^ to show its behavior essentially without 
finite-size effects. 

Finally, the spin- wave results for the von Neumann and 
Renyi entropies of the spin-1/2 antiferromagnetic Heisen- 
berg model on the square lattice are presented in Fig. |^. 
Fits to the form 5(i) = aL + b\nL + c give oi ~ 0.384 
for 5i = 5 and 02 ~ 0.191 for 52 after extrapolating 
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FIG. 4: (color online). Extrapolated values for a,b,c in fits 
of the entanglement entropies of Fig. ^ to the form S{L) = 
aL + 6 In L + c, indicated by dashed lines. Inset shows the 
scaling of a with the order of the Renyi entropy a. Solid line 
is a fit to aa = aooe"/" with Ooo ^ 0.095, w ~ 1.40. 



to account for the difference in fits for different numbers 
of points used. The coefficients 6 and c converge more 
slowly with the number of points used, but by taking the 
converged value of Oq, and using Sa{L) — UaL in a linear 
fit versus Ini, we find the extrapolated values shown in 
Fig. ^. Interestingly, it appears that the pref actor of the 
logarithmic term InL is essentially the same for all of 
the Renyi entropies including the von Neumann entropy, 
while the prefactor of the area-law term aL is a function 
of a; numerically, we find that the coefficients aa for the 



/ G 



inset 



a- Renyi entropies very closely follow Oq = 
of Fig. 1^. For the area-law coefficient of the single-copy 
entanglement entropy 5oo, this predicts aoo — 0.095. Ac- 
cording to spin-wave theory, therefore, the leading-order 
relation Si/Soo ^ 2 found in Rcf. |35| for critical spin 
chains does not appear to hold in two dimensions. We 
note that because of the additive logarithmic correction 
the entropies do not immediately saturate to a strict area 
law as is suggested by simply using small values of the 
staggered field h and extrapolating to h 0. An ad- 
ditional check on the validity of the spin-wave entangle- 
ment entropy is provided below using DMRG data for 
Heisenberg ladders. 



V. HEISENBERG LADDERS 

For a purely numerical approach to finding the entan- 
glement entropy of spin- wave Hamiltonians that is appli- 
cable to arbitrary couplings and boundary conditions, in 
particular to the Heisenberg ladder geometry of Ref. ||, 
we consider the Hamiltonian 



(25) 



where is the symmetric matrix describing the coupling 
between spins. Assuming the expansion about the clas- 



sical Neel state remains valid, the Dyson-Maleev trans- 
formation leads to the linear spin-wave Hamiltonian 



Hb = -BnocI + 



(26) 



with 



^Nccl = 



-(^V2)Ei,-/i. - hSN, 



{h + 

Bij = —SJij. We are thus left with the 
problem of diagonalizing the quadratic boson Hamilto- 
nian in Eq. (26). A and B are real symmetric matrices, 
and for our purposes we will assume that A ]s also diag- 
onal. The present approach is well-known,EjOL3tZl but 
here we describe the steps that lead to a rather straight- 
forward numerical implementation. Switching to first 
quantization, we use position and momentum variables 
in the Schrodinger representation through the standard 
relation b = {x + ip)/ \f2. Then the Hamiltoirian (|2^ ) 
becomes (written using vector notation) 

Hb = SnccI - ^Tr(A) + ^{p'^Gp + x'^Vx), (27) 

where G = A—B, V = A+B. We will assume that V and 
G are both positive-definite, but this must be checked for 
a given calculation. For the spin-wave Hamiltonian (|^) 
with nearest-neighbor coupling and PBCs, for example, 
the minimum eigenvalue is h for both G and V. We then 
seek the real, symmetric, and positive-definite matrix W 
such that WGW = V, in terms of which Eq. ( p7|) can be 
written as 



Hb = + -{p + iWxfG{p ~ iWx) . (28) 



Since G is positive-definite, if there is a state lipo) 
that satisfies {P — iWx)\ipo) = then j^/'o) is the 
ground-state of Eq. (Bq) with vacuum energy Eq = 
-EnocI ^ Tt{A — GWjj2. In the position basis we 
have {dxi + J2j^ij^j)i'o{x) = 0, and it can be eas- 
ily checked that the normalized solution is ipo{x) — 
[dct{W/TT)]^/^ cxp{—x'^Wx/2). The corresponding den- 
sity matrix po{x,x') — iPq(x)iPq{x')* is 



Pf){x,x') = [ dct 



W 



1/2 



exp 



- ]^{x'^Wx + x''^Wx') 



(29) 

The matrix W can be found by solving the generalized 
eigenvalue problem GVU = UVd in the form 



U^VU = Vd, UU'^ = G, 



(30) 



where Vd is a diagonal matrix. Then the inverse of W is 
given by 



(31) 



and satisfies W~^VW'^^ = G, so that W is found by 
inverting. If G is the identity as for the Klein-Gordon 
equation, for example, the conditions ( |30|) reduce to the 
usual problem of finding an orthonormal eigenbasis for 
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FIG. 5: (color online). Ground-state spin fluctuations J" 
(scaled by 4), von Neumann entanglement entropy S, and 
Renyi entropy S2 in the spin- 1/2 antiferromagnetic Heisen- 
berg model on 4A^ x A'^ ladders with OBCs. As shown in 
the inset, for each A'^ the subsystem is the left half, i.e., a 
2N X N region. Sohd (spin-wave) and dashed (DMRG) lines 
are guides for the eyes only. QMC results for T are shown as 
evidence for the convergence of the entropies, especially for 
N = 7,8. Cf Ref. | 



behaviors are similar. The subsystem sizes are 2N x N 
for each TV. The fluctuations are again close but slightly 
underestimated, while the von Neumann and Renyi en- 
tropies are over-estimated. Although it is difficult to 
compare numerical fits due to the small number of points 
for DMRG and even-odd oscillations, the DMRG data 
appear to approach the spin- wave results for larger N and 
fits of the spin- wave results to the form S — ai-f-folni-f-c 
give roughly a ~ 0.20 for the entanglement entropv and 
a ~ 0.10 for ^2, which are consistent with Fig. || after 
dividing by 2 for OBCs. The spin- wave theory does not 
perform as well for small N , and in particular does not 
adequateljy capture the gapped behavior of the even-A'^ 
ladders, £3 which is why there are much less oscillations. 
Again, the presence of the additive logarithmic correction 
implies that the entropies do not immediately saturate 
to a strict area law; contrary to previous interpretationsEl 
the DMRG data is also consistent with this fact, and 
QMC data for T shows that this is not due to lack of 
convergence. 



VI. CONCLUSIONS 



V , U"^ = U'^ , and W is simply the matrix square root 
of V. For the spin-wave Hamiltonian considered here, 
however, the more general procedure for finding W must 
be used. 

In the Schrodinger representation the number operator 
is given by hi = (xf + pf — l)/2, so that 

N ' 2 2N 2N ^ ' 

i 

As in Eq. (^3|) we impose the constraint e = S, but since 
in general the Hamiltonian is not translationally invari- 
ant the sublattice magnetization is zero only on average 
and may be non-zero on a given site. Application of 
Wick's theorem gives 

{h,h,) - {h,){h,) = + l{[W-% + (33) 

and the fluctuations can then be calculated using 
Eqs. (|T^), (pj|). Meanwhile, to compute the entangle- 
ment entropy of a region Q, let be the eigenvalues of 
XP where 

X^{x,x,}^^[W-%, P = {p,p^) = ^W,^ (34) 

restricted to Then the entanglement entropy is 

given by Eq. (p3[). 

Fig. [5] presents DMRG data and spin- wave results for 
47V X N Heisenberg ladders with OBCs, showing that the 



We have computed the von Neumann and Renyi en- 
tanglement entropies for the spin- 1/2 Heisenberg antifcr- 
romagnet on the square lattice using the modiflcd spin- 
wave theory for finite lattices. The results arc in good 
agreement with numerical results, in particular with the 
Renyi entropy 52 found from valence bond QMC calcula- 
tions and the entanglement entropy computed in DMRG 
for small ladders, and confirm the expected area law with 
important additive logarithmic corrections. We have also 
computed the spin fluctuations and found excellent agree- 
ment with QMC calculations, exhibiting an area law with 
a multiplicative logarithmic correction. These results 
support the view that linear spin-wave theory describes 
very well the ground-state of the spin-1/2 Heisenberg an- 
tiferromagnet on the square lattice. Interestingly, the es- 
sential features of entanglement entropy, namely the area 
law contribution and the important additive logarithmic 
correction, are captured at this level of approximation. 
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